Chaos Phenotypes in Fluids 

Y. Charles Li 



Abstract. I shall briefly survey the current status on more rigorous studies 
of chaos in fluids by focusing along the line of chaos phenotypes: sensitive 
dependence on initial data, and recurrence. 
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1. Introduction 

Chaos has two distinctive phenotypes: sensitive dependence on initial data, and 
recurrence. Sensitive dependence on initial data means that no matter how small 
the initial condition changes, after sufhciently long time, the change will reach order 
one. This phenotype can also be observed in non-chaotic systems, for example, in 
an explosive system. Together with the second phenotype, it can often identify 
chaos. Recurrence means that the orbit repeatedly re-visits the neighborhood of 
its initial point. 

Sensitive dependence on initial data can often be proved in a neighborhood of a 
homoclinic orbit or a heteroclinic cycle. The best technique for accomplishing this 
is the so-called shadowing lemma |14j . On the other hand, recurrence can often be 
established from general measure-theoretic arguments. 
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There are several good prototypes for understanding chaos in fluids: 2D Navier- 
Stokes equations under periodic boundary conditions and periodic forcing in both 
space and time, the Couette flow, the Poiseuille flow, and the plane Poiseuille flow, 
to name a few. When the so-called Reynolds number is large enough, all these 
flows develop chaos. When the Reynolds number just crosses its threshold, the 
chaos is often transient (that is, it has a flnite life time). With the increase of 
the Reynolds number, the life time of chaos increases too. Transient chaos is very 
poorly understood even for ordinary differential equations. The rigorously proved 
chaos is always eternal. Sometimes, one can observe intermittent chaos in fluids, 
that is, the chaos state and the regular state repeatedly exchange in time to eternity. 

In order to identify the chaos phenotypes in fluids, one can try to prove a re- 
currence theorem for Euler equations |17j [18j . and establish the existence of a 
heteroclinic cycle As a chaos phenotype, the recurrence theorem proved in 

[17j is quite satisfactory, while the heteroclinic cycle search in [11] is still at a 
conjecture stage. Another relevant topic to the dynamical system study of fluids 
is invariant manifold. It turns out that an unstable manifold can be proved for 2D 
Navier-Stokes equations, while invariant manifolds for Euler equations are a tough 
open problem [15j . All these studies were conducted on the simplest prototype — 
2D Navier-Stokes equations under periodic boundary conditions. For other pro- 
totypes mentioned above, very little is understood about chaos. Their analytical 
studies were often focused upon linear instability. There is a famous Sommerfeld 
paradox saying that Couette flow is linearly stable for all Reynolds numbers as flrst 
calculated by Sommerfeld [25j , but experiments show that Couette flow is unstable 
when the Reynolds number is large enough. This paradox lies at the heart of un- 
derstanding turbulence inside the inflnite dimensional phase space. An analytical 
study on the paradox is conducted in [20] . Another study was given in [9]. A 
complete resolution of the paradox will require a complete understanding of the 
dynamics at least in a certain region of the inflnite dimensional phase space. On 
the other hand, investigating the paradox may lead the way and have significant 
impact in the studies of fluid dynamics in phase spaces. Recently, there has been a 
renaissance in numerical dynamical system studies on the paradox [23] |22) [7] |29| 
[§] [26j [3] [27] d]. New non-wandering solutions (fixed points, periodic orbits, and 
relative periodic orbits) are discovered. Better estimates on the critical Reynolds 
number and perturbation threshold are obtained. It seems that pipe Poiseuille flow, 
plane Couette flow, and plane Poiseuille flow share some universal coherent struc- 
ture of "streak, roll, and critical layer" . The coherent structure seems quite robust 
in the sense that they survive for increasing Reynolds number. There seems to be 
also good agreement with experiments [5 . In the inflnite Reynolds number limit, 
the "streak" part of the flxed points tends to a limit shear (one of the inflnitely 
many shear flxed points of 3D Euler equations) ; and the "roll" and "critical layer" 
disappear [30] |28| . The limit shear is characterized by a condition |19| . 

The article is organized as follows. In section 2, we shall focus on recurrence. 
In sections 3-7, we shall discuss the issues related to sensitive dependence on initial 
data, especially the heteroclinics conjecture, linear spectra and their zero-viscosity 
limits, and the Melnikov integral. In section 8, we shall briefly discuss the Som- 
merfeld paradox. Section 9 is a discussion on future directions and open problems. 
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2. Recurrence 

The classical Poincare recurrence theorem can be stated as follows: 

Theorem 2.1 (Poincare Recurrence Theorem). Let fX, be a finite mea- 

sure space and f : X t-^ X be a measure-preserving transformation. For any 
E Cz T, (a-algebra of subsets of X), the measure 

^{{x &E\3N, f"{x) (^Eyn> N})^0 . 

That is, almost every point returns infinitely often. 

The geometric intuition of the Poincare recurrence theorem is that in a finite 
measure space (or invariant subset), the images of a positive measure set under 
a measure-preserving map will have no room left but intersect the original set 
repeatedly. The measure of the space X being finite is crucial. For example, 
consider the two-dimensional Hamiltonian system of the pendulum 

(2.1) i = y , y — — sinx . 

Its phase plane diagram is shown in Figure [TJ If the invariant region includes orbits 
outside the cat's eyes, then the measure of the region will not be finite, and the 
Poincare recurrence theorem will not hold. One can see clearly that the orbits 
outside the cat's eyes will drift to infinity. 




Figure 1 . The phase plane diagram of the pendulum equation. 

The investigation of recurrence of fluid dynamics in an infinite dimensional 
phase space encounters a serious problem that natural finite dimensional measures 
(e.g. Gibbs measure) do not have good counterparts in infinite dimensions. It 
is well-known that the kinetic energy and enstrophy are invariant under the 2D 
Euler fiow. But it is difficult to use them to define finite measures in infinite 
dimensions. It seems possible to study the Poincare recurrence problem directly 
from Banach norms rather than measures. Along this direction, the concept of 
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measure-preserving transformation, which is so critical in finite dimensions, is not 
necessary anymore. 

Consider the 2D Euler equation 

dtu + {u ■ V)m = ~Vp, V • M = 0; 

where u = (wi, U2) is the vefocity of the fluid, and p is the pressure. The boundary 
conditions are crucial in determining the structure of the dynamics in the phase 
space. The simplest boundary condition is the periodic boundary condition, that is, 
topologically the fluid motion is confined on the 2-torus . Under such a boundary 
condition, certain recurrence can be established as discussed below 17 . The 2D 
Euler equation is globally well-posed in H^{T'^) (s > 2). We also require that 

(2.2) / u = 0. 

The theorem stated below says that any H^{T'^) (s > 2) solution to the 2D Euler 
equation returns repeatedly to an arbitrarily small i7"(T^) neighborhood. 

Theorem 2.2. [17] For any u e H''{T'^) (s > 2), any S > 0, and any T > 0; 
there is a u* ^ H'^iT'^) such that 

F'^^^{u)eBl{u*) = {ueH\T^) I \\u-u*\\H0ij.)<5} 

where {wj} is an infinite sequence of positive integers, and is the evolution 
operator of the 2D Euler equation. 

Periodic boundary condition makes fluid dynamics a more natural dynamical 
system. The hope is that certain conclusions derived from the periodic boundary 
condition may be universal among all boundary conditions, or at least in the region 
away from boundaries. Of course, there are always special properties that are 
unique to the specific boundary conditions, for example, the example in [18j |21j of 
non-return near initial point is mainly due to the special slip boundary conditions 
on the boundaries of an annulus. 

3. Heteroclinics Conjecture 

2D Euler equation is an infinite dimensional Hamiltonian system [l]. Under 
periodic boundary condition, its Hamiltonian structure can be written in terms of 
Fourier series |12) . In a Hamiltonian system, heteroclinics often comes from the so- 
called separatrix. Such a separatrix represents the instability in the Hamiltonian 
system. The general scenario of chaos generation is that when the Hamiltonian 
system is under proper perturbations, chaos is created near the separatrix. Proper 
perturbations often represent a good balance between energy dissipation and energy 
input. For instances, going from 2D Euler equation to 2D Navier-Stokes equation, 
proper perturbations can be obtained by proper external forcings, e.g. spatially 
and temporally periodic forcing [llj . 

For the heteroclinics search, it is more convenient to write the 2D Euler equation 
in the vorticity form, 

(3.1) dtn + {^,n} = o, 

where is the vorticity which is a real scalar-valued function of three variables t 
and X = (xi, X2), the bracket { , } is defined as 

{/,5} = (a.J)(a,,5)-a./)(9.,5) , 
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where ^' is the stream function given by, 

the relation between vorticity and stream function 4" is, 

^ = dxiU2 - dx^ui = A* , 
We pose the periodic boundary condition 

fl{t, xi + 2tt, X2) — ^(t, xi, X2) — ^{t, xi,X2 + 27r/a), 
where a is a positive constant, i.e. the 2D Euler equation is defined on the 2-torus 

We propose the foUowing conjecture 

• The HeterocHnics Conjecture: In the Sobolev space H^{T'^) {£ > 3), for 
any fixed point Q of the 2D Euler flow having an unstable eigenvalue, 
there is a pair of heteroclinic cycles asymptotic to the two fixed points f2 
and —fl. 

The nature of "a pair of heteroclinic cycles" is motivated from the following 
symmetries: 

(1) n{t,Xi,X2) > Q.{t,-Xi,-X2), 

(2) n{t,xi,x2) — > -n{-t,Xi,X2), 

(3) ^{t,Xi,X2) > -fl{t, -Xi,X2), or Q{t,Xi,X2) > -H{t,Xi,-X2), 

(4) n{t,xi,x2) ^n{t,xi + ei,x2 + 02), ydi,02- 

The first symmetry allows us to work in an invariant subspace in which all the 
uJk's are real- valued. This corresponds to the cosine transform in (j3.2p . Take il = 
r cos xi as an example, the second symmetry maps the unstable manifold (assuming 
its existence) of the fixed point Fcosa;! into the stable manifold (assuming its 
existence) of — F cos xi . The third symmetry maps the unstable manifold of F cos xi 
into the unstable manifold of — Fcosxi. By choosing 61 ~ tt, the fourth symmetry 
maps the unstable manifold of Fcosxi into the unstable manifold of — Fcosxi. To 
maintain the cosine transform, the 9i and 62 in the fourth symmetry can only be 
TT and n/a. If there is a heteroclinic orbit asymptotic to Fcosxi and — Fcosxi as 
t — !■ —00 and +00, then there may be two corresponding to the unstable eigenvector 
and its negative. In fact, both may lie on certain sphere in the phase space due to 
the constraint by the invariants. Then the third symmetry generates another pair 
of heteroclinic orbit asymptotic to — Fcosxi and Fcosxi as i — s- —00 and +00. 
Together they form a pair of heteroclinic cycles. 
Using the Fourier series 

(3.2) J2 Wfce*('=i^i+"'=^^^) , 

fcez2/{o} 

where LU-k = one gets the kinetic form of the 2D Euler equation 

UJk= A{m, n) UJniUJn , 



k—m+n 



where 



1 



n\ + (0712)^ m\ + (0:7712)^ 
Denote by S the hyperplane 

Yj — {bj \ ujk — , V even k2} 



mi ni 

m2 712 
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We have the following theorem. 

Theorem 3.1. 11 Assume that the fixed point Q — rcosxi has a 1- dimensional 
local unstable manifold , and W" n S 7^ 0; then the heteroclinics conjecture is 
true, i.e. there is a pair heteroclinic cycles to the 2D Euler equation that connects 
f2 = Fcosa;! and — f2. 

Notice that the existence of invariant manifolds around the fixed point f2 = 
Fcosxi is an open problem. 

4. Numerical Verification of the Heteroclinics Conjecture 

Besides the symmetries mentioned in last section, we will also make use of the 
conserved quantities: kinetic energy E = ^ |fc|~^w^ (where |A;p = kf + a^fc|) and 
enstrophy S = '^uj'^, which will survive as conserved quantities for any symmetric 
Galerkin truncation, to help us to track the heteroclinic orbit. We will only consider 
the case that all the w^'s are real- valued (i.e. cos-transform). 

We make a Galerkin truncation by keeping modes: {\ki\ < 2, \k2\ < 2}, which 
results in a 12 dimensional system. We choose a = 0.7. After careful consideration 
of the above mentioned symmetries and conserved quantities {E = S ^ 1), we 
discover the following initial condition that best tracks the heteroclinic orbit: 

(4.1) W(o,i) = 0.603624 , = = 0.357832 , 

^{2,1) — ^(-2,1) = 0.435632 . 

We used fourth-order Runge-Kutta scheme. We also tested even higher-order 
Runge-Kutta schemes which do not improve the accuracy too much. Starting from 
this initial condition, we calculate the solution in both forward and backward time 
for the same duration of T = 11.8, and we discover the approximate heteroclinic 
orbit asymptotic to 2cosa::i and — 2cosxi as i — *■ —00 and +00, as shown in Fig- 
ure [21 Then the third symmetry generates another heteroclinic orbit asymptotic 
to — 2cosa;i and 2cosa;i as t ^ —00 and -foo. Together they form a heteroclinic 
cycle. Finally the second symmetry generates another heteroclinic cycle. That is, 
we have a pair of heteroclinic cycles. Notice also that the approximate heteroclinic 
orbit in Figure [2] has an extra loop before landing near — 2cosa::i. This is due to 
the fc2 = 2 modes in the Galerkin truncation. For smaller Galerkin truncations, 
the heteroclinic orbits can be calculated exactly by hand and have no such extra 
loop [13] [16] . and existence of chaos generated by the heteroclinic orbit can be 
rigorously proved in some case [16] . 

Remark 4.1. We have also conducted numerical experiments on Galerkin trun- 
cations by keeping more modes: {|fci| < 4, |fc2| < 4} and < 8, |/e2| < 8}. We 
found orbits that have similar behavior as the approximate heteroclinic orbit in 
Figure [51 but their approximations to heteroclinics are not as good as the one in 
FigurelH We also tested the heteroclinics conjecture for models of 2D Navier-Stokes 
equation |llj . 

5. Linear Instability and Invariant Manifold 

The 2D Navier-Stokes equation in the vorticity form can be written as 

(5.1) dtn + {^,n}^e[An + fit,x)], 
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Figure 2. The approximate heteroclinic orbit projected onto the 
(wj-i 0), ^(i,i))-plane in the case of the {|fci| < 2, \k2\ < 2} Galerkin 
truncation of the 2D Euler equation. 



where e = 1/R is the inverse of the Reynolds number, we will consider the same 
boundary condition as for (|3.ip . f{t,x) is the forcing, 

fdx = 0. 

T2 

For the external force / = Fcosxi, Q — Fcosxi is a shear fixed point, where 
r is an arbitrary real nonzero constant. Choose a G (0.5,0.84). There is a > 
such that when e > e*, the fixed point has no eigenvalue with positive real part, 
and when e G [0,e*), the fixed point has a unique positive eigenvalue [15] . Notice 
that this unique eigenvalue persists even for linear Euler (e = 0). In fact, for linear 
Euler (e = 0), there is a pair of eigenvalues, and the other one is the negative of 
the above eigenvalue. Precise statements on such results are given in the theorem 
below. Using the Fourier series p.2p where we work in the subspace where all the 
ojk's are real- valued, we get the spectral equation of the linearized 2D Navier-Stokes 
operator at the fixed point fl = 2 cos a; i. 



(5.2) 



where k e ZV{0}, p ^ (1,0), 



k-\-np ' 



A{q, r) 



1 



Ar^ = A{p, k 
1 



np) , and 





91 


ri 




92 


r2 



^ + (ar2)2 91 +(0:92)^. 

(In fact, the ^„'s should be counted twice due to switching q and r, but the differ- 
ence is only a simple scaling of e and A.) Thus the 2D linear NS decouples according 
to lines labeled by k. The following detailed theorem on the spectrum of the 2D 
linear NS at the fixed point ft = 2 cosxi was proved in 
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Theorem 5.1 (The Spectral Theorem [15]). The spectra of the 2D linear NS 
operator h5.2\) have the following properties. 

(1) (afc2)^ + (fci +"-)^ > I7 Vn G Z/{0}. When € > 0, there is no eigenvalue of 
non-negative real part. When e = 0, the entire spectrum is the continuous 
spectrum 

-ia\k2\, ia\k2\ 

(2) k2 — 0, ki — 1. The spectrum consists of the eigenvalues 

A = -en^ , n e Z/{0} . 

The eigenf unctions are the Fourier modes 

w„pe"^i + c.c. , yConp eC , n e Z/{0} . 

As e the eigenvalues are dense on the negative half of the real axis 

(— cx), 0]. Setting e = 0, the only eigenvalue is X — of infinite multiplicity 
with the same eigenf unctions as above. 

(3) k2 — —1, ki =0. (a), e > 0. For any a G (0.5,0.95), there is a unique 
£*("). 



V32 - 3a6 - 17a4 - 16a^ 1/1 



2(a2 + l)(a2 +4) ' (a^ + 1) 

where the term under the square root on the left is positive for a G 
(0.5,0.95), and the left term is always less than the right term. When 
e > e*(a), there is no eigenvalue of non-negative real part. When e = 
e*{a), X = is an eigenvalue, and all the rest eigenvalues have negative 
real parts. When e < e*(a), there is a unique positive eigenvalue A(e) > 0, 
and all the rest eigenvalues have negative real parts. e~"'^A(e) is a strictly 
monotonically decreasing function of e. When a G (0.5,0.8469), we have 
the estimate 



/a2(l-a2) a4(a2+3) 



2(a2 + l) 4(a2 + i)(a2 + 4) 



e(a2 + 1) < A(e) 



02(1-02) 

<V 2(a^ + l) ^ ' 

where the term under the square root on the left is positive for a G 
(0.5,0.8469). 



/02(l-02) o4(o2+3) ^ ^.^ ^^^^^ /o2(l 



y 2(a2 + l) 4(a2 + i)(a2 ^4) - ,_o+ y 2(o2 + 1) 

In particular, as t ^ 0^, A(e) = 0(1). 

(h). e = 0. When o G (0.5,0.8469), we have only two eigenvalues Aq 
and — Aq, where Aq is positive. 



y 2(a2 + l) " 4(o2 + i)(a2+4) ^ ° ^ y 2(a2 + 1) 
T/ie rest of the spectrum is a continuous spectrum [— io, ia\ 
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(c). For any fixed a e (0.5,0.8469), 
(5.4) lim A(e) = Aq . 

(4) Finally, when e — 0, the union of all the above pieces of continuous spectra 
is the imaginary axis iR. 

Based upon the above spectral theorem, the following invariant manifold the- 
orem can be proved. 

Theorem 5.2 (Invariant Manifold Theorem [E]). For any a e (0.5, 0.95), and 
e e (0, e*(a)) where 6,(0;) > satisfies iS.S]} . in a neighborhood of fl = 2cosxi in 
the Sobolev space iJ^(T^) (£ > 3), there are an 1- dimensional C°° unstable manifold 
and an l-codimensional C°° stable manifold. 

6. Zero- Viscosity Limit of the Linear Spectra 

One of the goals of the work |15j is to study the zero viscosity limit of the 
invariant manifolds of the 2D NS. For this study, it is crucial to understand the 
deformation of the linear spectra as e ^ O"*". Of course, studying this limit is of 
great interest in its own right. It is an interesting but difficult analysis problem 
too. We have conducted some numerical studies. We truncate (|5.2[) at different sizes 
and compute the eigenvalues of the resulting matrices. We increase the truncation 
size until we see reliability of the result. We also tested the continued fraction 
approach [12| for computing eigenvalues, the result is much worse. So we dropped 
the continued fraction approach. 

We find that the spectra of the linear Navier-Stokes and Euler operators can 
be classified into four categories in the zero viscosity limit 11]: 

(1) Persistence: These are the eigenvalues that persist and approach to the 
eigenvalues of the corresponding linear Euler operator when the viscosity 
approaches zero. (e.g. at 2D and 3D shears, and cat's eye.) 

(2) Condensation: These are the eigenvalues that approach and form a con- 
tinuous spectrum for the corresponding linear Euler operator when the 
viscosity approaches zero. (e.g. at 2D and 3D shears, cat's eye, and ABC 
flow.) 

(3) Singularity: These are the eigenvalues that approach to a set that is not in 
the spectrum of the corresponding linear Euler operator when the viscosity 
approaches zero. (e.g. at 2D and 3D shears.) 

(4) Addition: This is a subset of the spectrum of the linear Euler operator, 
which has no overlap with the zero viscosity limit set of the spectrum of 
the linear NS operator, (e.g. cat's eye.) 

When fci = and k2 — 1, a — 0.7, the unique e* in (|5.3p belongs to the interval 
0.332 < e* < 0.339, such that when e < e*, a positive eigenvalue appears. We test 
this criterion numerically and find that it is very sharp even when the truncation of 
the linear system (|5.2|) is as low as |n| < 100. As e 0+, we tested the truncation 
of the linear system (|5.2p up to |n| < 1024 for a = 0.7, the patterns are all the 
same. Below we present the case \n\ < 200 for which the pattern is more clear [llj . 
Figure [3fa) shows the case e — 0.14 where there is one positive eigenvalue and all 
the rest eigenvalues are negative. Figure[3^b) shows the case e = 0.13 where a pair 
of eigenvalues jumps off the real axis and becomes a complex conjugate pair. Figure 
(He) shows the case e = 0.07 where another pair of eigenvalues jumps off the real axis 
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(a) e = 0.14 



iOOO -5000 -4000 -3000 -2000 -1000 1000 
(b) e = 0.13 



0.4r 

0.3- 
0.2- 
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-0.3- 



-2500 -2000 -1500 -1000 -500 500 

Re(^) 

(c) e = 0.07 



1400 -1200 -1000 



-600 -400 -200 200 



(d) e = 0.03 



Figure 3. The eigenvalues of the hnear system (15.21) when fci = 
and ^2 = 1, a = 0.7, and various e. 



and becomes a complex conjugate pair. Figure [Sljd) shows the case e = 0.03 where 
another pair of eigenvalues jumps off the real axis and becomes a complex conjugate 
pair, while the former two pairs getting closer to each other. Figure ID^a) shows the 
case e = 0.0004 where many pairs of eigenvalues have jumped off the real axis and 
a bubble is formed. Figure [DJb) shows the case e = 0.00013 where the bubble has 
expanded. Including many other case testings, our conclusion is that: As e ^ 0"*", 
the limiting picture is shown in Figure [3I^c). Setting e = 0, the spectrum of the 
line fci = and fc2 = 1 of the linear Euler operator has been established rigorously 
(Theorem l5.ll) and is shown in Figurelljd), where the segment on the imaginary axis 
is the continuous spectrum. Comparing Figures lUJc) andHfd), we see that the two 
eigenvalues represent "persistence" , the vertical segment represents "condensation" , 
and the two horizontal segments represent "singularity" . Next we study one more 
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(a) e = 0.0004 



(c) e — > 0'^ limiting picture 
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-0.8 
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(b) e = 0.00013 
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Figure 4. The eigenvalues of the hnear system (15.21) when fci = 
and k2 — 1, a — 0.7, and various e (continued). 



hne: fci = and k2 = 2 {a — 0.7). In this case, there is no unstable eigenvalue. 
Figure O^a) shows the case e — 1.5 where all the eigenvalues are negative. As e is 
decreased, the eigenvalues go through the same process of jumping off the real axis 
and developing a bubble. Figure [5l^b) shows the case e = 0.00025 where the bubble 
has expanded. As e — > O"*", the limiting picture is similar to Figure |4jc) except that 
there is no persistent eigenvalue. The cases fci = and k2 > 2 {a — 0.7) are all the 
same with the case fci = and ^2 = 2 (a = 0.7). Figure [5l^a) shows the limiting 
picture of the entire spectrum of the linear NS operator as e ^ 0+ . Figure Elb) 
shows the entire spectrum of the linear Euler operator (e = 0) given by Theorem 

The fascinating deformation of the spectra as e — s- 0+ and the limiting spectral 
picture clearly depict the nature of singular limit of the spectra as e — > 0+ . In the 
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1.5r 



-5 -4 -3 
(a) e = 1.5 



■ -7 -6 -5 -4 -3 -2 

(b) e = 0.00025 



Figure 5 . The eigenvalues of the hnear system 
and ^2 = 2, a = 0.7, and various e. 



when ki — 



(a) e — > 0"^ limiting picture 



(b) € = 



Figure 6. The entire spectrum of the hnear NS operator (|5.2p 
when a = 0.7, e — > 0"*" or e = 0. 



"singularity" part of the limit, there is a discrete set of values for the imaginary 
parts of the eigenvalues, which represent decaying oscillations with a discrete set of 
frequencies. Overall, the "singularity" part represents the temporally irreversible 
nature of the e ^ 0+ limit, in contrast to the reversible nature of the linear Euler 
equation (e = 0). Many other fixed points were also studied in 
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7. Melnikov Integral 

We propose to use the constant of motion of the 2D Euler equation, 

G = [ Vl^dx- I \u\^dx 

to build a Melnikov integral for the corresponding 2D Navier-Stokes equation (|5.ip . 
The goal is to use the Melnikov integral as a measure of chaos, for example, around 
the line of fixed points ft = Fcosxi parametrized by T. G is a linear combination 
of the kinetic energy and the enstrophy. The gradient of G in is given by 

which is zero along the line of fixed points Q = Fcosa;! (a basic requirement for 
the convergence of the Melnikov integral). We define the Melnikov integral for the 
2D NS dSU) as 

M = ^ r f VnG[An + f{t,x)]dxdt 

/ [n + A-^n][An + f{t,x)]dxdt 

-oo Jt^ 

If the heteroclinics conjecture is true, then one can evaluate fl in the Melnikov 
integral along the heteroclinic cycle. 

As an example, we choose the external force 

(7.2) / = asini cos(xi + ax2) . 
Then the Melnikov integral (|7.ip has the expression 

(7.3) M = Ml + a^M^ + Ml sin(to + 0) , 
where 

. „ Ms . M2 

sm = 



Ml = ^ I I I {^ + A"il])A17 dxidx2dt , 



y/Ml+Ml y/Mf+Mi 



— 00 

/ / {il + A^^fl) costcos(xi + ax2) dxidx2dt 

-00 Jo Jo 

/+00 / a 1^2'JT 

/ / (Q + A~^n) sintcos(xi + 0x2) dxidx2dt . 

^00 Jo Jo 



For a numerical simulation, we evaluate f2(i) along the approximate heteroclinic 
orbit in Figure [5] with £7(0) given by (|4.ip . The time integral is in fact over the 
interval [—11.8,11.8] rather than (— (X),oo), which already gives satisfactory accu- 
racy. This is because that VG decays very fast along the approximate heteroclinic 
orbit in both forward and backward time. Direct numerical computation gives that 

Ml = -29.0977 , M2 = -0.06754695 , M3 = . 

Setting M = in (|7.3[) . we obtain that 

430.77741 



sin(io + 7r) = 
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Thus, when 

(7.4) \a\ > 430.77741 , 

there are sohitions to M = 0. The goal is to test whether or not (|7.4p is a criterion 
for the intersection of certain center- unstable and center-stable manifolds. In [llj . 
many such tests have been conducted. 



8. Sommerfeld Paradox 

The Sommerfeld paradox was originally about Couette flow, in fact, many other 
fluid flows exhibit the same paradox. The Couette flow is between two infinite 
horizontal planes (Figure [7]) where the upper plane moves with unit velocity. The 
dynamics of such a fluid flow is governed by the Navier-Stokes (NS) equations 
together with a unique boundary condition. Specifically 

(8.1) Ui^t + UjUi^j = + eu.jj , = ; 

defined in the spatial domain Doo = K x [0, 1] x E, where Ui {i = 1, 2, 3) are the 
velocity components, p is the pressure, and e is the inverse of the Reynolds number 
e — 1/R. The following boundary condition identifies the specific flow 

Ui{t,Xi,0,X3) = 0, Ui(t,Xi, 1,X3) = 1, 

(8.2) Ui{t,xi,0,X3) ^ Ui{t,xi,l,X3) ^ 0, (i = 2, 3). 




Figure 7. Couette flow. 



The classical Couette flow is given by 

(8.3) 1*1=2:2, 1*2 ii3 = 0. 

In the study [20j . we will focus not only on the neighborhood of the Couette flow 
(|8.3p . but also on the neighborhood of the following sequence of oscillatory shears 

(8.4) ui = U{x2) = X2 + — sin(4n7rx2), ( ^-r- < c < ) , U2 = W3 = 0. 

n \ 2 47r 47r / 

The sequence of oscillatory shears has some remarkable properties: 
lim ||[/(a:2) - a:;2||Lp = 0, lim \\U{x2) - X2\\w^:P ^ 0- 

n — >oo n — >oo 

That is, in the norms, the oscillatory shears can be regarded as perturbations 
of the Couette shear. One can regard the oscillatory shears as high frequency 
Fourier monomodes. In experiments, these oscillatory shears can be regarded as 
high frequency noises. On the other hand, in the W-^'P norms (i.e. vorticity's 
norms), these oscillatory shears are not perturbations of the Couette shear. 
A more precise statement of the Sommerfeld paradox is as follows: 

• Mathematically, the Couette shear is linearly and nonlinear stable for all 
Reynolds number R, in fact, all the eigenvalues of the Orr-Sommerfeld 
operator satisfy the bound A < —C/R |24j . 
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• Experimentally, for any R > 360 (where i? = ^ in our setting there 
exists a threshold amplitude of perturbations, of order 0{R~^) where 
1 < M < X depends on the type of the perturbations [10) . which leads to 
turbulence. 

A mathematically more comfortable re-statement of this experimental claim is as 
follows: For any fixed amplitude threshold of perturbations, when R is sufficiently 
large, turbulence occurs. For any fixed i?, when the amplitude of perturbations 
is sufficiently large, turbulence occurs. There may even be an asymptotic relation 
between such amplitude threshold and R. 

The main idea of our study is as follows |20j : The oscillatory shears (|8.4|) are 
perturbations of the Couette shear. As n ^ oo, they approach the Couette shear 
in L°° norm. They are linearly unstable under the 2D Euler dynamics [20j . This 
leads to the existence of an unstable eigenvalue of the corresponding viscous Orr- 
Sommerfeld operator, when the Reynolds number R is sufficiently large j20j . Notice 
that these oscillatory shears are not fixed points anymore under the Navier-Stokes 
dynamics. Nevertheless, they only drift very slowly. The important fact is that 
here the unstable eigenvalue of the Orr-Sommerfeld operator is order 0(1) with 
respect to e = 1/i? as e ^ 0+. Intuitively, this fact should lead to relatively long 
finite time nonlinear growth near the oscillatory shears (and the Couette shear). 
Such a growth will manifest itself in experiments as transient turbulence. Here 
the amplitude of the perturbation from the Couette shear will be measured by the 
deviation of the oscillatory shears from the Couette shear and the perturbation on 
top of the oscillatory shears. 

9. Future Directions and Open Problems 

In terms of analytical studies on chaos in fluids, the promising future direction 
will be the periodic boundary condition model where analysis can benefit from 
Fourier analysis. The heteroclinics conjecture is a challenging open problem on top 
of the open problem on the existence of invariant manifolds for 2D Euler equation. 
On the other hand, Galerkin truncations and numerical simulations all indicate that 
such a heteroclinics should exist. The current numerical capability is still limited in 
dealing such simulations. Of course, future improvement of numerical capability is 
always a hope. The question of how useful Melnikov integrals can be for predicting 
chaos in fluids, is still unclear. Once again, the flrst obstacle is the open problem on 
the existence of invariant manifolds for 2D Euler equation. 2D Euler equation may 
not have any invariant manifold at all. In principle, Melnikov integrals serve as good 
predicting tools for chaos only when the phase space structures are very special 
in the sense that the Melnikov prediction on the intersection of certain center- 
unstable and center-stable manifolds leads to the existence of certain homoclinics 
or heteroclinics. Here for the 2D Navier-Stokes equation, the phase space structures 
are not clear at all, for example, what are the fixed points or limit cycles? how 
close they are from the simple fixed point we are studying? 

The phase space structures of Couette flow, Poiseuille flow, and plane Poiseuille 
flow etc. also drew a lot of attentions recently. Most of the works are numerical. Due 
to the limitation of the numerical capability, analytical results are very precious. 
Analytically, the norms of the phase space play a fundamental role which is often 
missing from numerical simulations. It may be possible to have a better picture of a 
small region in the phase space by feedback studies between analysis and numerics. 
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It may also be possible to prove the existence of certain nontrivial fixed points or 
limit cycles. In finite dimensional dissipative systems, chaos is often associated 
with certain homoclinics or heteroclinics. It is not clear whether or not this is still 
true for fluid flows with boundary layers. There are some positive indications from 
recent numerical simulations [6] . 
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